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Abstract 

A new value b = 0.0205ib0.0017 fm^ of the slope of the neutron form factor GEnio^) 
at = compatible with deuteron properties has been extracted by using a linear 
relation between h and ^^(l + ff') we found for a class of nonlocal potential models 
having the experimental values of both r/5 and Q. Another model dependent value 
^MHKZ _ 0.0206 lb 0.0014 fm^, which is also compatible with deuteron properties, 
has been determined by applying " constrained" unitary transformations to the local 
MHKZ potential model. The sensitivity of a small changes in the experimental values 
used for ro and Q on the value obtained for h is also investigated. 

PACS: 13.40.Fn; 25.30.Bf; 21.10.Ft; 03.65.Nk. 



1 Introduction 



Measurements of the elastic scattering of electrons by deuterons lead in a direct 
way to the deuteron form factors as functions of the momentum transfer q^. 
These form factors are of interest in connection with the study of the structure 
of the deuteron and the determination of the electromagnetic form factors of 
the neutron [1]. 

Correlations between the slope h of the neutron electric form factor GEniff) at 
= and the deuteron root mean square (rms) matter radius ru for different 
NN-potential models have been used to determine r^) [2,3,4]. Berard et al. [2] 
found a hnear relation with a positive slope between h and for a class of 
potential models. They used their experimental data for the ratio R{q^) of 
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the deuteron to the proton electric form factors in the range of momentum 
transfers 0.05 < < 0.50 fm~^ and the experimental value b = 0.0189±0.0004 
fm^ given by Krohn and Ringo [5] to extract the value r£, = 1.9635±0.0045 fm. 
Allen et al. [3] found a linear relation with a positive slope between h and r|, for 
the radial wave functions for a single local potential model and its family of the 
wave functions produced by unitary transformations. The experimental data 
of -R(g^) of Berard et al. [2] and the experimental value b = 0.0199±0.0003 fm^ 
given by Koester et al. [6] have been used to determine the model dependent 
value To = 1.952±0.004 fm and the model independent one r£, = 1.948±0.023 
fm indirectly from the plotted straight line. They used also the experimental 
data of R{q^) of Akimov et al. [1] in the range 0.05 < < 0.50 fm~^ to 
determine td = 2.005 ±0.118 fm. In an analogous procedure Mustafa et al. [4] 
found a linear relation between b and rj^ for the radial wave functions for the 
MHKZ local potential model [4] and its family of the wave functions produced 
by unitary transformations. They used also the experimental data of R{q^) of 
Berard et al. [2] and the experimental value b = 0.0199 ± 0.0003 fm^ given by 
Koester et al. [6] to determine ro = 1.9546 ± 0.0021 fm. 

In this work, correlations between b and deuteron properties, i.e. the asymp- 
totic S-state amplitude As, the asymptotic D/S ratio rj, the rms radius ru, 
the quadrupole moment Q and the binding energy Ej,, have been used to de- 
termine b. Two linear relations between b and A|(l + ?7^) have been found 
for deuteron potential models, one with a positive slope and the other with a 
negative one. 

At large distances, outside the nuclear potential range, the ^Si- and the ^Di- 
state radial wave functions {u and w, respectively) behave asymptotically as 



where 7^ = —^E^ [8] and m is the reduced neutron-proton mass. If the 
proton and the neutron are at a distance r apart, the deuteron rms radius td 
may be written as [7] 
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The deuteron quadrupole moment Q may be also given by [7] 
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The correlation with v4|(l + 77^) has been used because the term (n^ + w'^) 
which has been involved in the definition of rjj proportional to + if 
we use the asymptotic forms of the radial wave functions u and w of Eq. (1) 
and ignore the centrifugal term [1 + 3/(7r) + 3/(7^r^)]. 

The aim of this work is to extract a new value of b which is compatible with 
deuteron properties. This value is extracted indirectly by using the linear 
correlation between b and v4|(l + 77^) of nonlocal potential models (having 
the experimental values of both and Q) produced by applying unitary 
transformations to standard nonrelativistic potential models. 



In Sect. 2 we will outline the method which has been used to calculate a model 
value of the slope b of the neutron form factor GEn{(f) at = 0. The new 
determinations of this parameter using the local and nonlocal potential models 
will be distinguished in Sect. 3. In Sect. 4 a simple method to calculate the 
parameter b using one single potential model will be obtained. Finally, the 
sensitivity of a small changes in the experimental values used for r^) and Q on 
the extracted value of 6, i.e. and -ttt, will be given in Sect. 5. 



2 Calculating model values of b 



In order to extract the slope b of the neutron form factor GEn{(f') at = 
from elastic electron-deuteron scattering, we start from the expression of the 
deuteron electric form factor in the nonrelativistic impulse approximation, 
neglecting two-body contributions [7] 
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where GEp{q^) is the proton electric form factor, GE^q^) and Gqi^q^) are electric 
charge and quadrupole form factors, respectively, which are given in terms of 
the nonrelativistic S- and D-wave functions by [7] 
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and 



/ 2 \ 

CQ{q') = 2ll^uw-^j J2(f)rfr, (6) 

where jo and j2 are the spherical Bessel functions of order zero and two, 
respectively. The factor r = {q^/Arrip) is called the Darwin- Foldy correction to 
the nucleon form factors [9], where nip = 4.5749098 fm~^ is the proton mass. 
Note that, Eq. (4) is usually used because Cq <^ Ce for small g^. 

The value of b for a particular potential model is insensitive to the choice of the 
proton electric form factor GEp{<f)- In this work we take the parameterization 
of Simon et al. [10] for in fm~^ 



2N 0.312 1.312 0.709 0.085 
GEp{q ) - 7" + T. 7^ ~ T. 7^ + T. 7^ ■ ( ' ) 
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For the evaluation of the charge form factor Csiq^) given in Eq. (5) we split 
the integration according to 



R 

CE{q') = / + Jo (f) dr + ACEiq') , (8) 



where AC£;(g^) is the asymptotic analytic correction to CE{q^) given by 



ACE{q') = j[u' + w')jo{'^-^]dr. (9) 

R 



The first integration of Eq. (8) is calculated numerically from r = to r = i?, 
where is the hard-core radius - if any - and R is chosen to be IG-t-r^ fm 
which is beyond the range of the nuclear potential. The integral in Eq. (9) is 
carried out analytically from r = i? to r = oo by using the asymptotic forms 
of the radial wave functions u and w of Eq. (1). The analytic formula of the 
ACEiq"^) is given in Ref. [4], for which the values of C^(q'^) are correct in the 
momentum transferred range < < 25 fm~^. In an analogous procedure 
the deuteron rms radius given in Eq. (2) and the deuteron quadrupole 
moment Q given in Eq. (3) are evaluated. 

Given CE{q^) for a particular potential model, then the experimental values 
of the ratio R{q^) of the deuteron form factor GEoiq^) to the proton form 
factor GEp{q^) can be used to determine the neutron form factor [2] from 
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and hence, the slope h of GEnil^) at = 0, 
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The experimental values of R{q^) of Simon et al. [11] are used in this work. 
This set of data covers the range of values 0.044 < < 4 fm~^. 

Various parameterizations have been used for GEn{q^) [2,12,13,14,15,16,17,18]. 
In order to calculate b in this work we parameterize GEn{q^) by a polynomial 
of order three in as 



which is agrees well for 0.044 < < 4 fm~^. The model dependency of 6, 
here, is a result of the model dependence of the charge form factor CE{q^) via 
the radial wave functions. 



3 The new determinations of b 

In order to extract a new value of the slope b of the neutron form factor 
GEn{q^) at = 0, a linear relation between b and v4|(l + rf) with a positive 
slope has been found for standard nonrelativistic potential models as shown 
in Fig. 1. These models and the "names" given to them are the potentials 

of Glendenning and Kramer "GKl, , GK9" [19], Lacombe et al. "PARIS" 

[20], Mustafa et al. "MHKZ" [4], Reid "RHC, RSC, RSCA" [21], Machledit 
et al. "MACH-A, -B, -C" [22], Machledit et al. "Bonn-F, -Q" [23], de Tourreil 
and Sprung "TS-A, -B, -C" [24], de Tourreil et al. "TRS" [25], Hamada and 
Johnston "HJ" [26], Mustafa and Zahran "MZ" [27], Mustafa "A, B" [28], 

Mustafa et al. "rl, r3, r7" [29], Mustafa "LI, L2, 1, 2, , 6" [33] and 

Mustafa "a, b, , i" [30]. The value 

6 = 0.0296 ± 0.0118 fm^ (13) 

which is compatible with the experimental values of As, rj and Eb but not 
with the experimental values of ro and Q is extracted from this correlation. 
It is the value of b corresponding to the experimental value A^{1 + rj^) = 
0.7817 fm~^. This value of b is greater than the previously extracted value 



GEniq'') = bq^ + cq^ + dq' 
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b = 0.0199 ± 0.0003 fm^ by Koester et al. [6] and its error Ab is relatively 
large. Therefore, the result obtained for b using the local potential models is 
not very reliable. 

Moreover, it is desirable to obtain a value of b which is not only compatible 
with the experimental values of As, rj and Eb, but also with the experimen- 
tal values of Q and ro, which have dependencies on the interior parts of the 
radial wave functions. The experimental values of these quantities are listed 
in Table 1 . Contributions to the value oivD and Q from meson exchange cur- 
rents and other relativistic corrections Ar^ = 0.0034 fm and AQ = 0.0063 
fm^ of Kohno [35] are neglected in the present work. Kermode et al. [36] and 
Mustafa and Hassan [37] have proven that Q and ro are not purely asymp- 
totic quantities by showing that phase-equivalent potentials can have different 
quadrupole moments and radii. Therefore, the linear relation between b and 
+ ?7^) obtained in Fig. 1 has been re-drawn using only potential models 
which reproduce the experimental values of both Q and r£,. These phase- 
equivalent potential models are produced by using "constrained" short-range 
unitary transformations like those used previously by Kermode et al. [36] for 
the radial wave functions {ui, Wi) of a local potential models. The transformed 
wave functions {ui,Wi) of a nonlocal potential model are given by 

oo 

Ui = Ui-2g{r) j g{s)ui{s)ds , 



oo 

Wi = Wi-2g{r) j g{s)wi{s)ds . (14) 



Various parameterizations have been used for the function g{s) [37,38,39,40,41]. 
We have chosen the parameterization g{s) = Cs{l — /5s)e~"'*, where s = 
r — Tf. and C = [4a;^/(a^ — 3/5a -|- SP"^)]^^"^ is a normalizing factor such that 
g'^{r)dr = 1. The parameters a and P are adjustable parameters. The pa- 
rameter a can be regarded as the "range" and (3 as the "strength" of the 
unitary transformation. The values of a < 1.5 fm~^ produce wave functions of 
different asymptotic behaviour [37]. In this work, the parameter a is assumed 
to changed from a = 1.5 fm~^ to a = 7 fm~^ in steps of 0.0001 fm~^. For 
each value of a, the parameter j3 is changed from (3 = fm~^ to (3 = 4 fm~^ 
in steps of 0.0001 fm~^. The function g{r) is of short range, hence, the wave 
functions u and w of a local potential model and u and w of the corresponding 
nonlocal potential model are the same in the asymptotic region. This implies 
that the asymptotic quantities As and Ajj and hence, r] = Ajj/As, are the 
same for both the local and nonlocal potential models, but td and Q could 
be different. 

The values of the non-locality parameters a and P are adjusted to produce 
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transformed wave functions having the experimental values of r^) and Q. This 
could not be achieved for all the local potentials considered. It is always pos- 
sible to find for a given value of a a value of (3 which corresponds either to 
the experimental value of r^) or to the experimental value of Q. In fact, for a 
given a one often finds two solutions for (3 as is shown by the two branches 
in Fig. 2. But it was not always possible for some potential models to find 
pairs such that both ru and Q could be reproduced. In this case we fix 

the non-locality parameters which give the experimental value of r^) by 

searching for the closest value of Q with respect to the experimental value or 
vice versa. If both values would be fitted, the two type of curves possess an 
intersection point as is demonstrated in Fig. 2 for selected potential models. 
The transformed wave functions having the experimental values of both 
and Q are compared to the radial wave functions in Fig. 3 for selected po- 
tential models. Altogether, we have studied forty-nine potential models. Out 
of these, we only found seventeen transformed models giving the experimen- 
tal values for both ro and Q. The pairs of the non-locality parameters 
having the experimental values of both ro and Q are listed in Table 2. 

In this case of applying "constrained" unitary transformations, a linear rela- 
tion with a negative slope between h and A\{l+7f) is found as shown in Fig. 4. 
The points representing the transformed potential models lie on or closely scat- 
tered around the straight line. We would like to mention that the five points 
representing the nonlocal potentials which are phase equivalent to the family 
of the Glendenning and Kramer potential models (GK2, GK3, GK5, GK7 and 
GK8) [19] lie on a separate line with a different negative slope (see also, Fig.l 
in Klarsfeld et al. [42] in which they used deuteron potential models in an 
empirical relation and found that the point representing the Glendenning and 
Kramer potentials also not lie on the empirical line, but lie on a separate line 
with a different slope). However, these old models predict a relatively larger 
deuteron binding energy Eb in agreement with the then experimental value 
Eh = —2.226 ± 0.003 MeV [43]. Therefore, these five potential models are not 
considered in Fig. 4. The new value 



5 = 0.0205 ± 0.0017 fm^ (15) 

is extracted from the straight line of Fig. 4 using only the remaining twelve 
nonlocal potential models. This value of b is compatible with the experi- 
mental values of As, rj, td, Q and E^. It is more consistent with the pre- 
viously published values b = 0.0189 ± 0.0004 fm^ of Krohn and Ringo [5] and, 
b = 0.0199 ± 0.0003 fm^ of Koester et al.[6]. The standard error in Eq. (15) 
is much less than that of Eq. (13), but is greater than those of the previous 
measurements. 
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4 Determination of b using MHKZ potential model 



So far, we have used two methods to extract the slope parameter b. The first 
uses the linear relation between b and of the local potential models 

as shown in Fig. 1. The second is based on Fig. 4 using the nonlocal potential 
models having the experimental values of both and Q. The two extracted 
values in Eqs.(13) and (15) differ substantially. 

A third method to extract b does not use the linear relation between b and 
yl|(l + ?7^) and a large number of potential models. It uses only one single 
potential model which is the potential model of Mustafa et al. [4]. The impor- 
tance of the MHKZ potential model comes from the fact that it reproduces the 
experimental values of As and t]. The unitary transformation applied to the 
MHKZ wave functions are constrained to produce a phase-equivalent trans- 
formed wave functions having the experimental values of both and Q. The 
model value, 



is extracted by using the phase-equivalent potential of the MHKZ potential 
model [4]. It is also compatible with the experimental values of As, rj, Eb, 
To and Q. This value is in very good agreement with the value obtained in 
Eq. (15) of this work. 



5 The sensitivity of a small changes in the experimental values 
used for td and Q on the extracted value of b 

The dependence of the determined value of b on the experimental values used 
for r£, and Q has been investigated. The whole procedure has been repeated 
twice, first by allowing for a small change 6rD in the value assumed as an 
experimental value of r^i and fixing the experimental values of As, rj, Eb and 
Q to find the corresponding change Sb. The published experimental values of 
To are in the range between r£, = 1.947 ± 0.029 fm of Akimov et al.[l] and 
rD = 1.9635 ± 0.0045 fm of Berard et al.[2]. The value of td of Akimov et 
al.[l], which is the lower limit of the experimental determination of ro, has 
been used in this work to study the correlation between the determined value 
of b, Eq. (15), and the value used for ro as an experimental value. The value 



,MHKZ 



0.0206 ± 0.0014 fm^ 



(16) 



6 = 0.0176 ±0.0013 fm^ 



(17) 
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is extracted in this case by using the transformed wave functions of potential 
models having the values of both r^i and Q. We show that the decrease in the 
value of leads to a corresponding decrease in the value of h. Therefore, we 
find 



: 0.6744 fm. 



5r 



D 



In the second case, only the experimental value used for Q is changed to 
obtain 6b/ SQ. The latest two published measurements of the experimental 
value of the deuteron quadrupole moment Q are the value of Reid and Vaida 
[44], Q = 0.2860±0.0015 fm^, and the value of Bishop and Cheung [31], Q = 
0.2859 ± 0.0003 fm^. To study the correlation between the extracted value of 
b and Q we used the experimental value Q = 0.2860 ± 0.0015 fm^ [44]. The 
result of the analysis is 



6 = 0.0205 ± 0.0011 fm^ (19) 

From these calculations we show that the increase in the value of Q leads to 
a small decrease in the value of b resulting in 



^ = -0.0052. (20) 
Similarly, we found 

-0.098 fm3/% (21) 



SAs 
and 



CI 

— = -0.002 fm^ . (22) 
or] 

From these results we conclude that the correlation between b and r^? is the 
dominant one. 

Now, we would like to point out the reason for the changing slope of the cor- 
relation between b and A|(l + rj'^) from being positive in the case of the local 
potential models of Fig. 1 to a negative value in the case of the nonlocal poten- 
tial models of Fig. 4. Since b is proportional to r^, and Td in turn proportional 
to v4|(l + ?7^), the value of b is proportional to A|(l + rj^). The values of As 
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and rj and then + 77^) are the same for the local and the correspond- 

ing nonlocal model, but b and r^) are different. We see that the value of b in 
the case of using the radial wave functions of the local model is considerably 
greater than that in the case of using the transformed wave functions of the 
corresponding nonlocal potential model. Therefore, the slope will change sign 
in the case of using the transformed wave functions of the nonlocal potential 
models. 
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Table 1 

The experimental values used in the analysis of this work. 



Value Ref. 

rz)=1.9547±0.0019 fm [9] 

Q= 0.2859±0.0003 fm^ [31] 

^5=0.8838±0.0004 fm-^/^ [32] 

r/=0.02713±0.00006 [33] 

Eb= -2.224575±0.000009 MeV [34] 



Table 2 

The pairs of the nonlocality parameters (q, /3) of the unitary transformations pro- 
ducing transformed wave functions having the experimental values of both ro and 
Q. 



local 


Ref. 


a 


P 


potentials 




fm-i 


fm~^ 


MHKZ 


[4] 


5.70890 


2.9672 


GK2 


[19] 


1.5818 


0.3966 


GK3 


[19] 


2.0516 


0.9938 


GK5 


[19] 


1.6421 


0.4143 


GK7 


[19] 


2.0455 


1.0127 


GK8 


[19] 


2.0894 


1.0066 


RSC 


[21] 


2.0650 


0.7787 


RHC 


[21] 


2.7131 


1.2650 


RSCA 


[21] 


1.8211 


0.7141 


MACH-C 


[22] 


1.5363 


0.6680 


rl 


[29] 


2.2160 


0.9669 


r3 


[29] 


2.3759 


1.1125 


r5 


[29] 


2.9344 


1.4908 


r6 


[29] 


3.5364 


1.8926 


r7 


[29] 


3.0106 


1.6228 


TRS 


[25] 


1.4821 


0.6472 


HJ 


[26] 


3.0968 


1.4556 
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Figure 1. The variation of the slope b versus ^^(l + ?7^) of standard local potential 
models. The middle (lower) part of the graph is magnified in the upper (lower) inner 
frame. The value of b of Eq. (13) is extracted from this straight line corresponding 
to A|(l + r/2) = 0.7817 fm-i. 
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1.6 2 2.4 2.8 3.2 3.6 4 1.6 2 2.4 2.8 



H (fm-^) fl[ (fm^) 

Figure 2. The variation of the "strength" parameter /? versus the "range" parameter 
a for some local potential models, indicated on the graphs, which have intersection 
points, i.e., have a pair (a, (3) which give the experimental values of both Q and r/j. 
The solid (dotted) lines represent the values of a and (3 which give the experimental 
values of Q (ru). 
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0.5 1 1.5 2 2.5 2 4 



r (fm) r (fm) 




r (fm) r (fm) 

Figure 3. Transformed wave functions having the experimental values of both ro 
and Q produced by the unitary transformation (solid line) are compared to the local 
potentials which are represented by the dotted line. The upper (lower) curves are 
the u (w) wave functions. 
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Figure 4. The correlation of the slope b versus + rj"^) of the twelve nonlocal 

"transformed" potential models having the experimental values of r/j and Q. The 
new value of b of Eq. (15) is extracted from this line. It is the value of b corresponding 
to the experimental value of ^1(1 + rj"^) = 0.7817 fm~^. 
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